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1 .  INTRODUCTION 


Atomic  oxygen  and  translational  temperature  are  involved 

in  important  chemical  and  collisional  processes  in  the 

thermosphere,  but  existing  methods  for  measuring  thermospheric 

profiles  of  the  0-atom  density  and  translational  temperature 

have  serious  deficiencies,  as  pointed  out  by  Zachor  and  Sharma 

(1989)  .  Sharm.a,  ct  .  (1987)  proposed  to  measure  the  two 

profiles  by  a  new  rem.ote  sensing  method  based  on  limb  radiance 

data  provided  by  a  spaceborne  Fabry-Perot  system  with  channels 

centered  at  147  am  and  63  urn.  wavelength,  corresponding  to  the 

intermultiplet  ground  electronic  state  01  transitions 
3  3  3  3 

{  Pq  and  report  is  the  second  of  two 

that  evaluate  the  feasibility  of  the  proposed  experiment 
concept.  The  earlier  one,  an  interim  report,  will  be  referred 
to  as  Part  1,  and  the  present  final  report  as  Part  2.  The 
reports  evaluate  two  fundamentally  different  approaches. 

Part  1  analyzed  a  single-channel  version  of  the  proposed 
technique,  in  which  the  limb  radiance  data  corresponds  to  just 
one  of  the  01  lines  spectrally  resolved.  It  v;as  tested  for 
both  the  147  umi  and  the  6  3  am  line,  using  synthetic  data 
corresponding  to  various  levels  of  signal-to-noise  and  spectral 
resolution.  The  147  am  line,  because  it  is  not  as  strongly 
absorbing  as  the  6  3  am.  lino,  proved  to  be  the  better  choice. 
Application  of  an  onion-peeling  technique  to  the  147  am 
synthetic  spectral  data  yielded  stable  temperature  and  0-atom 
density  solutions  from  300  km  down  to  an  altitude  between  130 
and  90  km,  depending  on  the  noise  level  and  spectral  resolution. 
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However,  it  was  determined  that  a  spaceborne  Fabry-Perot 
sensor  having  the  sensitivity  and  spectral  resolution  needed 
for  retrievals  down  to  90  km  altitude  would  require  rather 
large,  cooled  foreoptics  and  an  advanced  etalcn  design.  In 
particular,  the  wire  mesh  etalon  v;ould  have  to  be  much  larger 
in  size  than  we  believe  is  practicable  unless  the  Fabry-Perot 
system  is  confocal,  i.e.,  is  based  on  spherical  mesh  etalon 
reflectors.  While  a  confocal  etalon  providing  the  required 
high  resolving  power  (~5  x  10  )  is  theoretically’  possible, 
this  combination  of  design  and  performance  remains  to  be 
demonstrated . 

Here,  in  Part  2,  we  evaluate  an  alternative  approach  in 
which  the  data  consists  of  a  pair  of  limb  radiance  profiles 
corresponding  to  the  two  unresolved  01  lines.  This  approach, 
and  the  one  evaluated  in  Part  1,  will  be  referred  to  herein¬ 
after  as  the  nonresolved  and  resolved  approaches,  respectively. 
In  the  nonresolved  case,  recovery  of  the  translational  tempera¬ 
ture  is  based  physically  on  the  relative  01  line  intensities, 
that  is,  on  the  relative  upper-state  populations  of  the  two 
transitions  for  the  0-atoms  in  a  given  volume  element.  In  the 
resolved  case,  retrieval  of  temperature  relies  on  the  spectral 
shape  of  the  volume  emission,  that  is,  on  Doppler  broadening  of 
the  147  um  line. 

In  evaluating  the  nonresolved  approach,  we  inverted 
synthetic  limb  radiance  data  that  included  noise.  The  re¬ 
trievals  are  performed  by  an  algorithm  which  uses  the  onion¬ 
peeling  method  (see,  e.g.,  Goldman  and  Saunders,  1979)  over 


nost  of  the  300-to-90  km  altitude  range,  and  uses  the  global- 
fit  method  (see,  e.g.,  Carlotti,  1988)  near  200  km  altitude 
where  a  singularity  causes  divergence  in  the  onion-peeling 
method.  An  interesting  finding  of  the  study  is  that  changes 
in  the  01  limb  radiances,  typically  at  the  tangent  height 
200  km,  due  to  a  perturbation  in  0-atom  density  near  this 
altitude,  are  indistinguishable  from  changes  due  to  a  tempera¬ 
ture  perturbation;  i.e.,  the  four-elemient  Jacobian  matrix  of 
radiance  partial  derivatives  is  singular.  This  singularity, 
which  does  not  occur  at  other  tangent  heights,  is  more  easily 
avoided  in  the  global-fit  technique  than  in  the  onion-peeling 
method.  Based  on  the  retrieved  0-atom  density  and  temperature 
profiles  corresponding  to  various  levels  of  noise  in  the 
synthetic  data,  we  estir.ate  basic  system  requirements  for  a 
spaceborne  Fabry-Perot  sensor.  It  is  concluded  that  the  non- 
resolved  approach  is  cxper im.ental ly  feasible  using  currently 
available  technology. 
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2.  RADIATIVE  TRANSFER  EQUATIONS 


The  spectral  radiance  N^(z)  along  a  path  through  an 
emitting  and  absorbing  medium  with  boundary  condition 
Ny(0)  =  0  is  given  by 


z 


2 

J  n(x)  kv(x) 
y 


(y) 


dy/47l 


(1) 


where  n(z)  is  the  density  of  the  emitting/absorbing  species, 
k^{z)  is  the  absorption  coefficient,  and  ^(z)  is  the 
spectrally  resolved  volume  emission  rate,  defined  as  the 
number  of  photons  emitted  per  unit  time  per  unit  wavenumber 
interval  per  unit  volume.  The  exponential  factor  in  the 
integrand  is  the  spectral  transmittance  of  path  segment  yz; 
this  transmittance  times  ^(y)  dy/47t  is  the  contribution  of 
path  element  dy  at  point  y  to  the  spectral  radiance  at  point 
2.  By  applying  equation  (1)  to  a  path  through  the  limb  and 
evaluating  at  a  point  z  =  Z  above  the  atmosphere  (see 
Fig.  1),  one  obtains  the  limb  spectral  radiance.  In  general, 
ky  and  ^  depend  on  the  path  temperature  distribution  T{z)  and 
pressure  distribution  p(z).  We  assume  a  spherically 
stratified  atmosphere  characterized  by  the  vertical  profiles 
n(H)  (the  0-atom  number  density,  hereinafter  denoted  [Oj), 

T(H)  and  p(H).  A  limb  path  may  then  be  specified  by  its 
tangent  height  (Fig.  1).  A  spaceborne  sensor,  by  scanning 
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the  limb,  can  measure  the  limb  spectral  radiance  profile 
or  its  nonresolved  counterpart  N  (H^)  . 

In  the  present  work  N^(H^)  refers  to  the  spectral 
radiance  arising  from  a  single  01  transition  at  wavenumber  Vq 
(wavelength  Xq  =  147  |im  or  63  njn)  .  Associated  with  the 
transition  is  an  Einstein  coefficient  A,  upper  state 
statistical  weight  g^^  and  upper  state  energy  Values  of  A, 

g^^  and  for  the  two  01  lines  are  given  ^in  Table  1.  The 
corresponding  volume  emission  rate  assuming  LTE,  is  given 
by 

^  =  A.[0]  f  (v-VQ,T,p)g^  exp(-C2EyT) /Q{T)  (2) 

where  f(v, T,p)  is  the  normalized  Voigt  shape  of  the  line  due 
to  Doppler  and  collisional  broadening,  =  1.4388  K/cm  ^  is 
the  second  radiation  constant,  and  Q(T)  is  the  partition  sum, 

Q(T)  =  exp(-C2  X  226. 5/T)  +  3  exp(-C2  x  158. 5/T)  +  5  (3) 

The  absorption  coefficient  is  given  by 

k^  =  f(v-Vo,T,p)  S(Vo,T)  (4) 

where  S  is  the  line  intensity, 

3(V3,T)  =  [A/ {87CcVq^)  ]  [  (l-y) /ylg^  exp  (-C2E^/T) /Q  (T)  (5) 
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witn 


y  =  exp(-C2Vo/T) 


(6) 


Equation  (2)  for  the  volume  emission  rate  can  be  written  in  the 
alternative  form 


^  =  47t[0]ls,J,Q 


(7) 


where  the  source  function  J^g  is  2cVq^Y/ ( 1 -y)  .  Note  that  J^q 
is  equal  to  Planck's  function  for  the  photon  emission  rate,  at 
wavenumber  Vg,  of  a  blackbody  at  temperature  T. 

The  atmosphere  is  represented  as  a  series  of  concentric 
shells,  of  vertical  thickness  1  km,  within  which  [0],  T  and 
p  are  constant.  The  integrals  in  equation  (1)  can  be 
evaluated  exactly  for  the  contribution  of  each  limb  path 
segment  Az^  corresponding  to  a  given  shell  (Fig.  1) .  The 
resulting  expression  for  the  limb  spectral  radiance  is 


Nv(Ht)  = 


N 


N 


exp 


i  =  l 


1-exp  (-AXvi) 


JvOa 


(8) 


j=i  +  i 


where  =  [0]^  k^^  Az^  is  the  optical  thickness  of  path 

segment  Az^.  A  derivation  of  equation  (8)  from  basic 
principles  can  be  found  in  Bullitt,  et  al.  (1985) .  They  give 
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in  addition  an  expression  for  the  absorption  coefficient  that 
is  applicable  under  general  (LTE  or  non-LTE)  conditions. 

Ccllisicnal  broadening  has  a  very  small  effect  on  the 
integrated  01  line  radiances,  which  in  the  present  application 
correspond  to  limb  tangent  heights  >  90  km.  Even  so,  we 
used  the  Voigt  rat.her  than  Doppler  lineshape  for  f  (V)  . 

However,  in  calculating  the  collisional  (Lorentz)  width  we 
used  a  fixed  pressure  profile,  specifically,  the  one  in  the 
U.3.  Standard  Atmosphere  (1976).  With  this  procedure,  the 
147  |im  and  63  jim  limb  radiance  profiles  computed  from  the 
above  equations  depend  solely  on  the  vertical  profiles  [0]  vs. 
H  and  T (H) .  Retrieval  of  the  [0]  and  T  profiles  involves  a 
systematic  fitting  of  computed  147  ^m  and  63  fim  limb  radiances 
to  the  corresponding  measured  limb  radiances. 

It  should  be  noted  that  large  errors  in  the  retrieved  [0] 
and  T  profiles  would  result  if  one  were  to  assume  an 
optically  thin  limb  path,  which  would  correspond  to  evaluating 
equation  (8)  in  the  limit  as  the  Axs  approach  zero.  For 
example,  the  spectral  transmittance  at  the  center  of  the  63  |im 
line,  between  the  tangent  point  and  the  observer  at  z  =  Z,  is 
typically  0.007  for  tangent  height  =  120  km,  and  0.95  for 
tangent  height  =  240  km.  The  corresponding  transmittances 
at  the  center  of  the  147  |im  line  are  -^0.31  and  ~  0.98, 
respectively.  If  the  objective  or  the  radiative  transfer 
calculation  was  merely  to  estimate  or  N  ,  then  the 

error  resulting  from  the  optically  thin  assumption  would 
probably  be  judged  acceptable  above  =  240  km  and 
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unacceptable  below  this  tangent  height.  However,  even  the 
small  error  in  N  for  H.j,  »  240  km  would  cause  a  very  large 
error  in  the  retrieved  [0]  and  T  profiles,  as  we  shall 
demonstrate . 


3.  RETRIEVAL  METHODOLOGY 

The  retrieval  method  to  be  described  depends,  physically, 
on  differences  in  the  local  volume  emission  rate  for  the  two 
01  lines.  It  has  the  capability  to  transform  a  pair  of 
nonresolved  limb  radiance  profiles  N  for  the  two  lines 
into  a  pair  of  vertical  profiles  for  the  nonresolved  volume 
emission  rate, 

^(H)  =  A„*[0]g^„exp(-C2E^„/T)/Q(T);  n  =  1,2  (9) 

which  is  the  rate  given  by  equation  (2)  after  integrating  over 
the  spectral  line.  In  equation  (9),  n  =  1  corresponds  to  the 
147  Jim  line  and  n  =  2  indicates  the  63  jim  line;  [0]  and  T 
correspond  to  the  same  altitude  as  One  can  obtain  T  (H) 

from  the  ratio  ^1/^2  ‘ 

R  =  =  0.0653  expl-97,84/T  (H)  ]  (10) 

where  we  have  used  the  values  in  Table  1  for  the  A„,  g„  and 
^un"  using  the  inferred  T(H)  and  or  then 

determine  [0] .  Our  retrieval  method  does  not  solve  explicitly 
for  (H)  and  ^2  obtains  T  and  [0]  more  directly,  as 
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described  below.  The  method  depends  nonetheless  on  the 
sensitivity  of  the  ratio  R  to  temperature  for  the  recovery 
of  T{H),  and  on  the  sensitivity  of  and  ^2 

recovery  of  [0]  vs.  H.  The  value  97.84K  appearing  in  equation 
(10)  is  the  temperature-equivalent  difference  in  upper  state 
energies,  equal  to  Cj  • 

From  equation  (10)  we  find  that  the  error  in  temperature 
due  to  a  relative  error  dR/R,  if  T  were  determined  from  this 
equation,  is 

dT  =  (T^/97.84)  dR/R  (11) 

The  error  dR  would  originate  from  noise  in  the  limb  radiance 
data  and/or  inevitable  small  errors  in  the  solutions  for 
and  ^2*  Erom  equation  (11),  the  value  dR/R  =  0.01  results  in 
an  error  dT  =  4K  when  T  =  200K;  the  relative  temperature 
error  is  only  4/200  =  0.02.  If  dR/R  =  0.01  and  T  =  1200K,  a 
temperature  that  can  occur  near  300  )cm  altitude,  the  error  dT 
is  147K;  the  relative  error  is  0,12.  One  can  show  similarly 
that  the  relative  error  d[0]/[0]  is  roughly  three  times  dR/R 
or  d^/^.  Thus,  it  can  be  anticipated  that  temperatures  can  be 
recovered  more  accurately  in  the  lower  thermosphere  than  at, 
say  300  1cm,  where  the  translational  temperature  will  be  large 
compared  to  the  difference  (~  98K)  between  upper  state 
energies;  and  that  the  0-atom  density  can  be  recovered  to 
higher  relative  accuracy  than  the  temperature  over  most  of  the 
thermosphere.  Also,  it  is  clear  that  the  optically  thin 
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assumption,  while  it  would  cause  errors  in  R  less  than  a  few 
percent  above  240  km,  would  lead  to  rather  large  errors  in  the 
recovered  temperature  above  this  altitude. 

The  retrieval  algorithm  used  over  most  of  the  90-300  km 
altitude  range  is  based  on  the  onion-peeling  method,  which  is 
applied  one  layer  at  a  time  beginning  at  the  top  of  the 
atmosphere.  The  solution  values  [0]  and  T  for  the  layer  at 
altitude  H  =  are  defined  as  those  for  which  the  function 

2 

F([0],T)  =  ^  {  [Nn(HT)-Nn(HT)  ]/Nn(HT)  }2  (12) 

n=l 

is  a  minimum.  Again,  subscript  n  distinguishes  between  the 
147  |im  and  63  ^im  lines;  denotes  a  computed  limb 

radiance  (integral  of  the  spectral  radiance),  while  N^(H^)  is 
the  corresponding  measured  radiance.  Thus,  F([Oi,T)  is  a  sum- 
square  normalized  difference  between  the  predicted  and 
measured  limb  radiances  for  the  01  lines.  Of  course,  Nj^(H^) 
includes  the  net  effects  at  the  observation  point  z  =  Z  of  all 
layers  above  altitude  H  =  for  which  solutions  [0],T  have 
already  been  obtained.  The  procedure  used  to  start  the  onion¬ 
peeling  process  and  to  estimate  10]  and  T  for  layers  above 
the  starting  altitude  (330  km)  will  be  described  subsequently. 

The  pair  of  least-squares  normal  equations  obtained  by 
formally  minimizing  F([0],T)  are 
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2 

^3Nn(HT)  /9ak[Nn(HT)  -  Nn  (Ht)  ]  /Nn2  (%)  =  0;  (13) 

n=l 

k  =  1,  2;  tti  =  [0]  ,  02  =  T 

A  solution  [0],T  of  equations  (13)  is  found  using  Newton's 
method.  There  is  no  solution  when  the  determinant  of  the 
Jacobian  matrix  3N/3a  is  zero. 

The  term  "layer, "  as  used  in  the  above  description  of  the 
onion-peeling  procedure,  refers  to  specific  index  shells, 
one  )cm  in  vertical  thickness,  whose  altitudes  differ  by 
AH  >  1  km.  For  example,  the  tangent  height  in  equation 
(13) ,  which  is  also  the  altitude  of  an  index  shell,  will  take 
on  the  value  150  km  if  was  160  km  in  the  previous  cycle, 
since  Ail  is  defined  to  be  10  km  for  »  130  km.  A  constraint 
on  the  solutions  of  equation  (13)  is  that  T  and  log[0]  vary 
linearly  with  altitude  between  the  two  index  layers,  and 
connect  with  the  [0],T  solutions  for  the  higher  index  layer. 
Thus,  the  onion-peeling  cycle  in  which  is  150  km  actually 
yields  [0],T  solution  values  for  the  ten  shells  at  H  =  159, 

158,  ...,  150  km.  The  guess  values  required  by  Newton's 
method  are  based  on  the  presumption  that  the  slope  of  T(H)  and 
log[0]  vs.  H  will  not  change  from  one  cycle  to  the  next.  The 
retrieval  algorithm  uses  AH  =  4  km  between  130  and  110  km 
altitude,  and  AH  =  2  km  below  110  km.  These  values  are 
adequate  to  reproduce  the  structure  of  typical  [0]  and  T 
vertical  profiles.  Note  that  the  derivatives  3N^(H.j)/3a^  will 
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generally  become  smaller  as  AH  is  reduced,  which  will  increase 
the  potential  for  instability  in  the  retrieval  solutions  due 
to  noise  (and/or  propagation  of  noise  and  error) . 

The  recovered  temperature  profile  is  smoothed  above 
130  )cm  altitude  by  convolving  it  with  a  rectangular  window 
function.  Smoothing  is  performed,  in  each  onion-peeling 
cycle,  over  the  20-km  altitude  range  addressed  in  the  two  most 
recent  cycles.  Thus,  each  portion  of  the  final  T(H)  solution 
above  130  km  has  been  convolved  twice  with  the  boxcar,  and 
exhibits  no  slope  discontinuities. 

The  Jacobian  referred  to  above  was  evaluated  as  a 
function  of  tangent  height  for  two  model  [0],T  profiles. 

These  are  the  same  models  used  to  generate  synthetic  limb 
radiance  data  for  testing  the  retrieval  algorithm.  The 
finite-difference  Jacobian  elements  were  based  on  perturbed 
[0],T  profiles  that  differed  from  the  model  profiles  only 
between  H,j  and  H.p  +  lO  km,  and  that  represented,  over  this 
10-km  region,  linear  T  vs .  H  and  linear  log  [0]  vs.  H.  In 
other  words,  the  numerical  Jacobian  values  typified  those 
which  would  occur  in  applying  the  onion-peeling  algorithm  to 
the  synthetic  limb  radiance  data.  The  Jacobian  was  found  to 
be  singular  very  near  H.j,  =  200  km,  but  at  no  other  for 
both  models. 

It  is  well  known  that  the  onion-peeling  method  tends  to 
propagate  noise.  Hence,  the  T  and  [0]  profiles  recovered 
from  noisey  data  can  differ  substantially  from  the  true 
profiles  near  200  km  altitude,  and  consequently  the  Jacobian 
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singularity  can  occur  at  an  altitude  somewhat  higher  or  lower 
than  200  km.  This  uncertainty  in  the  location  of  the 
singularity  tended  to  defeat  the  simple  strategy  of  locating 
the  index  layers  well  away  from  the  200-km  level,  e.g.,  one  at 
220  km  with  the  next  lower  one  at  180  km.  That  is,  successful 
solutions  were  a  chance  occurrence.  However,  it  is  notable 
that  the  onion-peeling  procedure,  whenever  it  did  not  fail 
because  of  the  singularity,  produced  solutions  down  to  ~  90  km 
altitude.  In  fact,  the  errors  in  the  retrieved  T  and  [0] 
profiles  near  ~  100  km  altitude  were  generally  much  smaller 
than  those  above  ~  180  km.  It  is  apparent  that  the  effects  of 
noise  propagation  in  the  onion-peeling  method  are  offset  by 
(a)  an  increased  sensitivity  of  the  01  line  radiances  (ratio 
R)  to  temperature  at  the  lower  altitudes,  where  the  atmosphere 
is  cooler,  and  (b)  higher  signal-to-noise  (a  brighter  limb)  at 
the  lower  tangent  heights. 

As  an  alternative  to  the  onion-peeling  method,  we 
considered  the  approach  in  which  predicted  01  radiances  are 
fit  to  the  measured  radiances  simultaneously  for  several 
tangent  heights.  By  fitting  at  M  tangent  heights,  it  is 
possible  to  recover  [0]  and  T  simultaneously  for  N  ^  M 
altitudes.  This  method  of  retrieval  is  described  by  Carlotti 
(1988),  who  calls  it  the  global-fit  approach,  and  by  Chang  and 
Weinreb  (1985) ,  who  refer  to  it  as  the  Levenberg-Marquardt 
method  after  the  algorithm  they  used  to  effect  the  nonlinear 
least-squares  fit.  We  found  that  the  global-fit  method 
required  considerable  time  on  the  CDC  Cyber  860  computer  when 
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it  was  used  to  recover  [0]  and  T  for  just  a  few  levels  below 
the  starting  altitude;  e.g,,  600  sec  was  required  when 

M  =  N  =  3.  Indeed,  Carlotti  (1988)  describes  a  retrieval  of 
stratospheric  ozone  concentrations  by  the  global-fit  approach 
that  required  substantial  execution  time  on  a  Cray  computer. 

The  algorithm  finally  adopted  in  the  present  study 
performs  onion-peeling  above  280  km  altitude  and  below  180  km. 
The  global-fit  method  is  employed  for  the  retrieval  of  [0]  and 
T  between  these  two  altitudes.  Specifically,  we  use  a 
modified  Levenberg-Marquardt  algorithm  (subroutine  UNLSF  in 
the  International  Math  and  Statistics  Library)  to  fit 
simultaneously  the  measured  limb  radiances  at  =  230  and 
130  km  and  to  obtain  [0]  and  T  at  these  two  altitudes  (i.e., 
M  =  N  =  2),  with  the  constraint  that  log[0]  vs.  H  and  T (H)  are 
linear  between  280  and  230  km  and  between  230  and  180  km.  The 
function  that  is  minimized  is  +  F2,  where  F^  and  F2  are 
defined  as  F  (Eq.  12)  evaluated  for  =  230  km  and  for 

=  180  km,  respectively.  This  procedure  proved  effective  in 
avoiding  the  singularity  at  200  km  altitude  for  a  wide  range 
of  noise  levels  in  the  synthetic  data. 

The  starting  procedure  exploits  the  fact  that  above 
250  km  the  translational  temperature  is  nearly  constant,  while 
the  0-atom  density  tends  to  fall  off  exponentially  with 
altitude.  Assuming  these  properties,  and  also  that  the  limb 
above  =  250  km  is  optically  thin,  one  can  show  that  N  (H,p) 
is  approximately  exponential  with  the  same  scale  height  as  [0] 
vs.  H;  see  Sharma  et  al.  (1988) .  In  the  starting  procedure. 


a  single  radiance  scale  height  Hg  is  estimated  from  the  two 
limb  radiance  profiles  above  =  250  km,  and  is  used  to 
extrapolate  the  profiles  from  =  300  (the  largest  value  for 
which  synthetic  data  was  generated)  to  =  330  km.  The 
least-squares  normal  equations  (13)  are  then  solved  for  [0] 
and  T  at  H  =  330  km,  using  the  two  limb  radiances  for  this 
tangent  height.  The  constraint  on  the  solution  in  this 
initial  application  of  the  normal  equations  is  that  T  is 
constant  and  that  d  log[0]/dH  equals  the  constant  -1/Hg  above 
330  km  altitude.  Thus,  the  solution  defines  T (H)  and  [0]  vs. 

H  for  H  ^  330  km.  Note  that  the  optically  thin  assumption, 
while  it  is  used  in  this  starting  procedure  to  equate  the 
scale  heights  of  [0]  and  the  limb  radiance  data  N(H.j),  is  not 
invoked  in  computing  the  predicted  limb  radiances  N (H^)  in 
equation  (13) . 

Generally,  there  is  some  error  in  the  T  and  [0]  values 
obtained  at  H  =  330  km  due  to  the  assumptions  that  T  is 
constant  and  [0]  is  exponential  above  250  km.  The  subsequent 
onion-peeling  cycles  tend  to  produce  compensating  errors  that 
alternate  in  sign;  i.e.,  the  retrieved  profiles  exhibit 
oscillations.  The  oscillations  usually  decrease  with  each 
cycle  and  are  very  small  when  the  altitude  H  =  280  km  is 
reached.  This  is  the  reason  for  starting  the  retrieval  at 
330  km  and  for  choosing  280  km  as  the  highest  level  addressed 
by  the  global  fitting  procedure.  After  the  T  and  [0] 
profiles  are  recovered,  values  above  H  =  290  km  are  discarded. 
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In  summary,  the  retrieval  is  started  at  330  km  and  is 
based  on  the  described  onion-peeling  procedure,  with  a  spacing 
of  10  km  between  index  shells,  down  to  280  km.  The  global-fit 
technique,  with  index  levels  at  230  km  and  180  km,  is  used  to 
recover  T  and  [0]  between  280  and  180  km.  The  onion-peeling 
method  is  used  below  180  km,  with  AH  =  10  km  above  130  km;  AH 
is  4  km  between  130  and  110  km  and  is  2  km  below  110  km.  Each 
10-km  section  of  the  temperature  profile  above  130  km  is 
smoothed  twice  regardless  of  which  method  is  used  to  retrieve 
the  section. 


4.  INVERSIONS  OF  SYNTHETIC  DATA 

Two  models  for  the  vertical  profiles  of  0-atcm  density 
and  translational  temperature  over  the  altitude  range  90  to 
300  km  are  shown  in  Fig.  2.  These  models,  designated  ATOI2 
and  ATOI4,  were  used  to  compute  limb  radiance  profiles  for  the 
147  pjn  and  63  Jim  01  lines,  which  are  shown  in  Fig.  3.  The 
limb  radiance  profiles,  after  addition  of  Gaussian  noise, 
constitute  the  raw  synthetic  data  to  which  a  smoothing 
procedure  and  then  the  retrieval  algorithm  was  applied. 

The  sensor  system  that  would  measure  the  limb  radiance 
profiles  (see  Section  5)  is  envisioned  to  have  a  larger 
instantaneous  f ield-of-view  (IFOV)  above  160  km  than  below 
this  tangent  height  in  order  to  provide  a  system  noise- 
equivalent  radiance  (NER)  five  times  smaller  above 
H^  =  160  km.  In  other  words,  spatial  resolution  is  traded  for 
sensitivity  at  the  higher  tangent  heights,  where  N(H,p)  is  less 
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structured  and  is  failing  off  rapidly.  It  is  assumed  that  the 
instrument  samples  the  01  limb  radiances  at  intervals  of  2  km 
in  tangent  height.  Thus,  we  generated  synthetic  profiles  at 
2-km  intervals,  added  noise  samples  corresponding  to  a  given 
NER  value  (one-fifth  this  value  above  =  160  km),  and 
finally  smoothed  the  profiles  above  =  160  by  convolving 
them  with  the  11-point  cubic  polynomial  smoothing  function 
given  by  Savitzky  and  Golay  (1964) .  The  11  points  of  the 
smoothing  function  correspond  to  a  total  width  of  20  km.  The 
convolution  was  performed  five  times  in  succession,  starting 
each  time  at  a  different  tangent  height,  to  effect  smoothing 
above  160,  200,  220,  240  and  250  km,  respectively.  Smearing 
of  the  limb  radiance  profile  by  the  finite  instrument  IFOV  was 
not  accounted  for  per  se.  However,  above  160  km  the  smoothing 
procedure  produces  considerably  more  smearing  than  would 
result  from  the  system  IFOV.  Note  that  the  retrieval 
procedure  uses  only  the  limb  radiances  which  correspond  to  the 
i.ndex  levels  defined  earlier. 

Retrievals  were  performed  for  the  two  synthetic  data  sets 
(corresponding  to  the  models  ATOI2  and  AT0I4)  for  a  wide  range 
of  NER  values.  The  NER' s  used  for  the  63  )im  limb  radiances 
were  always  0.6  times  those  used  for  the  147  pm  limb 
radiances,  corresponding  to  the  relative  detector/preamp- 
limited  noise  levels  for  the  two  sensor  channels  (see 
Section  5) .  Three  retrievals  were  performed  for  each 
combination  of  atmospheric  model  and  NER  level;  these 
correspond  to  three  independent  sets  of  Gaussian  random 


-17- 


numbers  for  the  noise  samples.  We  will  present  results  only 
for  the  highest  NER  levels  that  resulted  in  stable  solutions, 
since  it  is  these  levels  that  define  the  required  sensor 
performance . 

Figure  4  shows  the  actual  and  retrieved  [0],T  profiles 
for  atmospheric  model  AT0I2  when  the  system  NER,  for 

^  160  km,  is  lO'^^W  cm'^sr”^  in  the  147  ym  channel  and 
6  X  10  ^^W  cm  ^sr  ^  in  the  63  ym  channel.  This  is  the  highest 
noise  level  for  which  stable  solutions  could  be  obtained.  In 
fact,  in  one  of  the  six  cases  (6=2  models  x  3  noise  sets) 
corresponding  to  this  noise  level,  the  solution  diverged  at 
160  km  altitude;  the  other  five  solutions  were  stable  down  to 
104  km  (the  altitude  of  the  peak  of  the  [0] -profile)  or  lower. 
Figure  5  shows  the  actual  and  retrieved  [0],T  profiles  for 
model  ATOI4  for  the  same  NERs. 

Figures  6  and  7  show  typical  results  for  models  AT0I2  and 
AT0I4,  respectively,  when  the  NER' s  are  half  the  above  values, 
i.e.,  5  X  10  and  3  x  10  cm  ^sr  respectively,  for  the 
147  ^im  and  63  |am  channels.  In  all  six  cases  corresponding  to 
this  noise  level,  the  solutions  were  stable  down  to  90  km 
altitude . 

The  four  solutions  display  similar  characteristics:  The 
errors  in  retrieved  temperature  tend  to  be  largest  at  the 
higher  altitudes,  specifically,  where  the  actual  temperature 
is  greater  than  600K.  Even  though  these  errors  and  noise  are 
propagated  by  the  onion-peeling  procedure,  the  temperature 
solutions  are  remarkably  stable  at  the  lower  altitudes.  In 
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contrast,  the  errors  in  the  retrieved  0-atom  density  profile 
are  significant  only  below  the  peak  of  the  [0]  profile  . 

These  errors  are  the  result  of  noise  and  noise/error 
propagation  in  combination  with  the  very  high  opacity  of  the 
63  |im  line  below  120  km  tangent  height. 

It  may  be  noted  that  high  signal-to-noise  levels  are 
required  for  successful  retrievals.  For  example,  the  NER 

“13  -2  “1 

value  of  5  X  10  W  cm  sr  for  the  147  channel,  which 
resulted  in  the  [0],T  solutions  shown  in  Fig.  7,  corresponds 
to  S/N  ^  400  at  tangent  height  160  km,  based  on  the  signal 
levels  given  in  Fig.  3.  The  NER  value  of  3  x  10  cm”^sr'^ 
at  63  [im  corresponds  to  S/N  2  x  10^  at  =  160  km.  The 
latter  S/N  value  is  probably  very  much  higher  than  needed  for 
successful  retrievals;  it  results  from  the  fact  that  the 
signal  (limb  radiance)  is  much  higher  at  63  Mm  than  147  ^im, 
while  the  detector-limited  noise  levels  that  can  be  achieved 
at  the  two  wavelengths  using  Ga:Ge  detectors  are  comparable. 

5 .  COHCLaSIONS 

Application  of  the  described  retrieval  method  to 
synthetic  OT  limb  radiance  data  has  established  that  vertical 
profiles  of  the  atomic  oxygen  density  and  translational 
temperature  can  be  successfully  recovered  over  the  altitude 
range  90  km  to  300  km  if  the  noise-equivalent  radiance  for  the 
147  Mm  channel  is  no  greater  than  5  x  10  cm'^sr”^  for 
tangent  heights  .4^  S  160  km,  and  no  greater  than  lO'^^W  cm~^ 
sr  ^  above  160  km.  It  is  assumed  that  the  63  M^i  channel  can 
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achieve  NERs  equal  to  0.6  times  these  values,  based  on  the 
fact  that  commercially  available  Ge:Ga  photoconductive 
detectors  can  provide  a  noise-equivalent  power  (NEP)  of 
3  X  10  at  63  )im  wavelength  while  a  stressed  detector 

of  the  same  type  can  provide  NEP  =  5  x  lO’^^W/H^"^^  at  147  |im. 

An  additional  requirement  is  that  the  sensor  resolve  the  limb 
radiance  profile  to  at  least  2.5  km  tangent  height  for 

<  160  km.  A  tangent  height  footprint  five  times  larger 
would  suffice  above  160  km  tangent  height. 

Design  goal  data  requirements,  i.e.,  noise-equivalent 
radiances  and  instantaneous  fields-of-view,  are  given  in 
Table  2.  The  listed  NERs  are  approximately  three  times  lower 
than  the  maximum  allowable  values  stated  above.  The  listed 
IFCVs  correspond  to  the  requisite  minimum  tangent  height 
resolutions  if  the  sensor  altitude  is  300  km.  The  lower  NERs 
for  H.J,  =*■  160  km  could  be  achieved,  for  example,  by  coaddi.ng 
the  outputs  of  detectors  in  a  five-by-five  mosaic. 

Figure  8  is  a  block  diagram  showing  the  principal 
components  of  a  two-channel,  non-scanning,  cryogenic  Fabry- 
Perot  sensor  to  measure  the  01  limb  radiance  profiles.  Design 
specifications  for  this  sensor  concept  are  listed  in  Table  3. 
The  values  given  will  provide,  according  to  our  calculations, 
the  design  goal  NER' s  in  Table  2.  Note  that  the  collecting 
telescope  is  20  cm  diameter  and  that  the  entire  system 
operates  at  a  temperature  of  lOK;  the  detectors  must  operate 
at  4K  or  lower.  Coaddition  of  the  25  detector  outputs  is 
assumed  for  tangent  heights  greater  than  160  km.  Also,  the 
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noise  bandwidth  of  1  Hz  applies  to  a  single  pair  (63  |jjn  and 
147  |jjn)  of  limb  radiance  samples.  This  implies  roughly  100 
seconds  for  measurement  of  the  two  limb  radiance  profiles, 
which  is  net  an  unreasonable  value.  Of  course,  tradeoffs  are 
available  among  the  parameters  of  Table  3,  which  represents 
a  baseline  model  for  engineering  design  of  the  sensor.  It  is 
also  possible  that  detector  NEP's  lower  than  those  in  Table  3 
will  be  available  in  the  future. 

The  baseline  model  represents  a  system  of  reasonable 
size,  with  components  that  are  well  wit.hin  current  state-of- 
the-art  technology.  It  is  concluded  that  the  experiment 
proposed  by  Sharma  et  al  (1937),  using  the  nonresolved 
measurement  approach,  and  using  the  retrieval  methodology 
described  herein,  is  technically  feasible. 
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Table  3.  Sensor  Baseline 
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FIGURE  CAPTIONS 


1. )  Coordinates  and  parameters  involved  in  calculating  the 

limb  spectral  radiance  (see  text) . 

2. )  Two  models  for  the  vertical  profiles  of  atomic  oxygen 

density  and  translational  temperature.  They  represent  a 
cool  thermosphere  (model  AT0I2,  represented  by  dashed 
cur’  3s)  and  a  relatively  warm  thermosphere  (model  AT0I4, 
full  lines)  . 

3. )  01  limb  radiance  profiles  calculated  for  the  two  models 

shown  in  Fig.  2.  The  dashed  lines  correspond  to  model 
AT0I2,  the  full  lines  to  model  AT0I4. 

4. )  Actual  (full  lines)  and  retrieved  (dashed  lines)  profiles 

of  [0],T  for  atmospheric  model  AT0I2.  The  noise- 
equivalent  radiances  for  the  147  pun  and  63  ^im  channels 
are  10’^^  W  cm'^sr”^  and  6  x  10’^^  W  cm'^sr  \ 
respectively . 

5. )  Same  as  Fig.  4,  except  the  atmospheric  model  is  AT0I4. 

6. )  Same  as  Fig.  4,  except  the  NERs  for  the  147  jim  and  63 

channels  are  5  x  10  W  cm  ^sr  ^  and  3  x  10 
“2  -1 

W  c.m  sr  ,  respectively. 
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7.)  Same  as  Fig.  6,  except  the  atmospheric  model  is  AT0I4. 


8.)  Sensor  block  diagram  for  the  nonresolved,  two-color 
concept . 
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FIGURE  8 
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